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ABSTRACT 

We investigate the relation between the star formation history and the evolution of the 
double-degenerate (DD) population in the thin disc of the Galaxy, which we assume to have 
formed 10 Gyr before the present. We introduce the use of star- formation contribution func- 
tions as a device for evaluating the birth rates, total number and merger rates of DDs. These 
contribution functions help to demonstrate the relation between star-formation history and the 
current DD population and, in particular, show how the numbers of different types of DD are 
sensitive to different epochs of star formation. 

Analysis of the contribution functions given by a quasi-exponentially decline in the star- 
formation rate shows that star formation from to 8 Gyr after thin-disc formation dominates 
the present rates and total numbers of Hen-He DDs and CO-i-He DDs. Similarly, the current 
numbers of CO-i-CO and ONeMg-i-X DDs come mainly from early star formation (< 6 Gyr) 
, although star formation from 4 to 8 Gyr continued to contribute more CO-i-CO than Hen-He 
DDs. The present birth-rates for CO-i-CO and ONeMg-i-X DDs are strongly governed by recent 
star formation (i.e. 8 — 9.95 Gyr). Star formation from < 7.5 Gyr does not contribute to the 
present birth rates of CO-i-CO and ONeMg-i-X DDs, but it has a distinct contribution to their 
merger rates. 

We have compared the impact of different star- formation models on the rates and numbers 
of DDs and on the rates of type la (SNIa) and core-collapse supernovae (ccSN). In addition 
to a quasi-exponential decline model, we considered an instantaneous (or initial starburst) 
model, a constant-rate model, and an enhanced-rate model. All were normalised to produce 
the present observed star density in the local thin disc. The evolution of the rates and numbers 
of both DDs and SNIa are different in all four models, but are most markedly different in the 
instantaneous star-formation model, which produces a much higher rate than the other three 
models in the past, primarily as a consequence of the normalisation. 

Predictions of the current SNIa rate range from w 2 to 5 x 10^^ yr^^ in the four models, 
and are slightly below the observed rate because we only consider the DD merger channel. 
The predicted ccSN rate ranges from 1.5 to 3 century"^, and is consistent with observations. 
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1 INTRODUCTION 

We know little about the star formation (SF) history of the Milky 
Way disc. Of the few ways available to explore the SF history, 
studying the stellar age distribution by the observation of a sample 
of different types of stars and comparing observed and synthetic 
colour-magnitude diagrams would be the most favored approach. 
The simplest view of the SF history would be a single star burst tak- 
ing place at the formation of the disc. Substantial evidence, how- 



* E-mail: syu@arm.ac.uk (SY) 
f E-mail: csj@arm.ac.uk (CSJ) 



ever, shows the disc has experienced a more complicated history, 
including repeated star bursts, an exponential declining SF or ex- 
tende d periods of enhanced SF ( lMaiewskilll993l : lRocha-Pinto et al.l 
1 20001) or some combination of all of these. 

The present SF rate is better understood. ISmith et alj J1978l) 
concludes the SF rate in the Galaxy is perhaps 5 M0yr~'^, of which 
74% take place in spiral arms, 13% in the interarm region, and 13% 
in the galactic center. Observations of Lyman continuum photons 
from O stars in giant H II regions give the SF rate in the di sc to be 
4.35 M 0yr~\ A SF rate Ri3.6 Afoyr"^ was suggested bv lMcKeel 
1 119891) from an analysis of thermal radio emission in HII regions 
around massive stars. The measurement of the mass of '^^Al in 
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the Galaxy implies a SF rate «4 Moyr"! foiehl et al.ll2006l) . The 
derivation of the SF rate is strongly associated with the initial mass 
function (namely the distribution of the mass of the new-born stars) 
and the birth rate of core collapse supernovae (type Ib/c and type 
II). 

Double degenerates (DDs) are a type of exotic binary consist- 
ing of two white dwarfs. As the evolved remnants of stars, the for- 
mation of close DDs requires their progenitors to have undergone a 
mass transfer stage, either common envelope or stable Roche lobe 
overflow phase. Ultra-compact DDs are especially interesting since 
not only can they be observed as active sources of electromag- 
netic radiation, e.jj. AM Canum Venaticorum star s, JCropper et al.l 
1 1 9981: iRamsav et ai]|2005l : iNelemans et"Zll2004h . but also the fi- 
nal products of DD evol ution may include t ype la supernovae 



JYungelson & Lividl2000l)). Theore tical r esults JEvans et alj|l987l : 
iNelemans et al.ll200ll : iRuiter et al.ll201(]| : IYu & Jeffervll20Iol) sug- 
gest that close DDs should also be significant sources of gravita- 
tional wave radiation, and be detectable by the space gravitational 
wave (GW) detector LISA (Laser Interferometer Space Antennae). 

The census of DDs therefore is important because they allow 
us to test the endpoints of stellar evolution. We need to know the 
present-day numbers, birth rates, merger rates and galactic distri- 
bution of DDs, and also the distribution of their properties {i.e. 
mass, orbital separation, chemical abundance and age). To inter- 
pret these, we need to investigate how the history of star formation 
in the Galaxy affects the DD population. We can achieve this using 
a population synthesis approach. 

In this paper, we study the correlation of the SF history and the 
evolution of DD population. It may help interpret the distribution 
of the physical properties of DDs as a function of the SF history. In 
fj2l we describe the approach to simulate a population of the DDs 
in the thin disc and discuss the details of how a quasi-exponential 
declining SF rate influences the birth rate, merger rate and total 
number of DDs. In ij3] we show the distributions of some important 
physical parameters of the DDs from different epochs of SF. In ^ 
we compare the impact of different SF models on the rates and 
numbers of DDs. We discuss the results for supernovae rates in 
different SF models in ij5] Some points are discussed in ^ and 
conclusions are drawn in ^ 



2 POPULATION SYNTHESIS OF THIN DISC DDS IN 
THE MONTE CARLO APPROACH 

In this section, we describe the population synthesis method to ob- 
tain a sample of DDs comparable with current observations, and the 
structure of the thin disc in which the DD population is distributed. 
The individual stella r-evolution tr a cks were co mputed using the 
method described bv iHurlev et al.l <200Ct l2002r) . Po pulation syn- 
thesis was carried out using the method described bv lYu & Jeffervl 
( 12010.) in an initial study of the present Galactic DD population. 



2.1 The star-formation contribution function 

We define a time representing the age of the thin disc to be t, such 
that the disc formed at t = 0, and the current age of the thin disc 
is tdisc- During the evolution of the thin disc, we assume that star 
formation takes place over an interval t = tst (0 ^ isf ^ idisc) 
to f = tat + Stat. We simulate the evolution of all binaries formed 
during this interval by computing a sample (denoted by *) of k 
primordial MS-l-MS (MS: main sequence) binaries with a total mass 
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Figure 1. The orbital periods (log Porb^day) of the DDs against the primary 
masses (Afp, left panel) in the MS-hMS binary progenitors and the time 
(tMS. right panel) which is from a primordial MS-hMS binaiy to when its 
descendant DD was just bom. 



MJj. This gives a mean sample SF rate S, = M^/5tsi. 5tsf is 
assumed to be small; thus all primordial MS-l-MS binaries formed 
during this interval evolve as though they formed at time igf. The 
binary-star evolution tracks are computed as a function of t over 
the interval log(t/yr) = 7.4 - 10.13 with 5 log t = 0.00867. 

We then compute the contribution of star formation during the 
interval St^i to the birth-rate, total number and merger-rate of DDs 
at time t — tdisc (the present epoch). This interval is denoted by 
the quantity A = tdisc — tst. 

For a given age of the thin disc tdisc, star formation in the thin 
disc runs from tsf = to tdisc- For convenience, we use 5 log tsf — 

5 log tdisc 

Two additional time variables are used to establish the stage of 
evolution of DDs in the sample. tMS represents the lifetime from 
birth of an individual binary system to the formation of a DD. tuo 
represents the lifetime from formation of a DD to its mergeiQ 

The birth-rate contribution is established by counting the num- 
ber of DDs which form during the interval tdisc to tdisc + ^tdisc- 
These are stars for which 



tdisc ^ tsf + tMS < idisc + <5tdisc. 



(1) 



Normalising by the mean sample SF rate 5** gives the fractional 
number 5n,,ncw(A) of DDs bom per unit mass of star's formed 
between t = tsf and tsf + Stsf. 

The merger-rate contribution is established in an analagous 
way. Counting stars for which 



idisc ^ isf + ^MS + ioD < tdisc + <5tdis 



(2) 



gives the fraction 5n*,mor(A). 

The total-number contribution 5n*,dd(A) represents the frac- 
tional number of DDs in the sample which exist at time tdisc- These 
stars are indentifed by 



tsf + ^MS < idisc ^ isf + tyiS + tDD- 



(3) 



^ We make the assumption that the global properties of the MS-hMS binary 
population from each star formation episode are the same; i.e each sample 
adopts the same initial parameter distribution. 

^ We have assumed the angular momentum evolution of a DD to be gov- 
erned only by gravitational wave radiation and mass transfer <Yu & Jeffervl 
I2OI0I) , although tidal intera ction and magnetic braking may play an impo r- 
tant role jMarsh et al.l2004HGokhale et alj2007l ; lFarmer & Roelofsl2010l) . 
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Table 1. The variables in the population synthesis simulation and their re- 
lationship. 



at time fdisc- The integral of the rate contribution function tells us 
the birth rate and the merger rate of DDs. 



time 




tdisc 


the age of the thin disc 


tsf 


the time of star formation (SF) 




(0 < ht < tdi.c) 


A = idisc - isf 


delay time 




from SF to the disc age 


Thin disc 




"(tdiac) 


total number 


"(tdisc) 


number rate 


C(A) 


contribution function to n 


C(A) 


rate contribution function to n 


5(trf) 


SF rate 


Sample^ 




5n,(A) 


number at tdisc 




from SF at t^f — >■ tsf + (5tsf 


a (A) 


contribution function 


a (A) 


rate contribution function 


S*(tsi) 


mean SF rate in the sample 


t All the values normalized by 5* (t^f ) . 





Identities 

C'.(A) = an4A)/ai,f 
Ri<5n.(A)/5tsf 
a(A) = |^^^ 

C'{A) = C'.{A)-5(ief) 
C'{A) = a(A)-5(tef) 

"(tdisc) = /o""'^ C-Cidisc - t)dt 
= /o*d- n(t)dt 

"(tdisc) = /o'"°'' (^(idisc - t)di 



These three quantities, let us call them star-formation num- 
ber contribution functions Srit, comprising (5n,,ncw, Srit^dd and 
5n*,mer, may then be combined with a model for the star formation 
history of the thin disc S{t) to establish the total current birth-rate, 
merger rate and number of DDs. 

From above, we have 5n,(A) to be the SF contribution at 
time isf to the number per unit SF rate at time isf , integrated over 
an interval 5tsf. Hence C* = dn^/dtsf ~ 5n,/Sts{ is the sample 
contribution function per unit SF rate per unit time at time f^f . 

To obtain total numbers in a real thin disc, we must multiply 
by the thin-disc SF rate to obtain the total contribution rates: 



C ^C- S{t, 



5n, 



S{U, 



(4) 



C also reflects the distribution of age of the DDs. For example, if 
the thin-disc SF rate is S{t), star formation at tst ~ will con- 
tribute (5n, (tdisc)/<5tsf ■ S(0) DDs with age tdisc; star formation at 
tsf ~ 1 Gyr will contribute (Jriast (tdisc — l)/i5fsf ■ ^(l) DDs with 
age tdisc — 1 Gyr, and so on. 

Since we know from Eq.|4]that star formation at time tsf con- 
tributes C (tdisc — tsf) DDs at thin disc age tdisc, we define a rate 
contribution function 

C = dC/dUisc (5) 

to be the contribution from star formation at tst to the number rate 



2.2 Birth rate, merger rate and total number of DDs 

We define n{t) to be the total number of DDs at thin disc age t, 
where n may represent new-born rincw, merged rimer or existing 
ridd DDs. Star formation from t^f to t^f + St^f will generate C ■ Staf 
DDs with age tdisc — tst. Hence, the overall number of new-born, 
alive or merged DDs at t = tdisc is 



(6) 



n(tdisc) = / ^(tdisc — tst)dtst 

Jo 

tdisc 

« 22 5n, {tdisc — tst) • S(tsi) 

tdiac 
^ ^ C'-.Jtsf. 



Then we are able to calculate the number rate ri (tdisc) of DDs 
from the contribution function C 



^l^t^discj — 



C(tdisc — tst)dtsf 



E 



(Sn, (tdisc —tst) 



<5ts 



S'(ts 



(V) 



£ C • 5ts 



where n represents the birth rate v or merger rate C,. We will give 
the details of these functions in the present simulation in '\2A\ 

Furthermore, since we know the overall rate of change in the 
number DDs i/(t) — C(t)> we can also calculate the total number of 
DDs in the thin disc at time tdisc (ridd (tdisc)) by the integral 



ndd (tdisc) = / \v{t) - C(t)]dt 
/o 



£[Kt)-C(t)]-'5t. 



(8) 



Consequently, we have two different methods to compute the 
current number of DDs in the thin disc. Eq.[6]represents the sum of 
contributions from each individual SF epoch to the total number by 
counting which DDs exist at the current epoch. Eq.[8]represents the 
integral of the birth-rate minus the merger-rate, or nett birth rate, 
over the entire SF history of the galaxy. The two methods should 
give the same result but, due to the limited grid in tst, the evalua- 
tion of the sum in Eq .[6] gives slightly higher numbers than Eq.[8] 
In the present simulation with quasi-exponential SF rate, the differ- 
ences are 1.57%, 0.98%, 0.51%, and 0.46% for He-l-He, CO-l-He, 
CO-l-CO, and ONeMg-l-X DDs, respectively. 

We list the time variables, numbers (rates), contribution func- 
tions and their relation in the thin disc and our simulation in Table[T] 



2.3 Input to the Monte-Carlo simulation 

The SF rate is assumed to be given by the quasi-exponential func- 
tion 



^(tst) = 7.92e-(*=f)''" +0.09(tst) Moyr"' 



(9) 
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Figure 2. The variation of the contribution of star formation tgf to the 
present birth rate and merger rate of DDs, C, in the quasi-exponential star 
formation rate model. In this figure, top panels represent the contribution to 
the present birth rate, and bottorn panels represent the contribution to the 
merger rate. See il2.4l for the details. 



where r = 9 Gyr JYu & Jeffervl uOlOf ). which prod uces « 3.5 
MQyi ~^ at the current epoch. This is consistent with IPiehl et al.l 
( l2006h . 

The contribution function for each starting epoch (tst) is ob- 
tained from a sample of primordial binaries. Each binary is defined 
by at least five initial parameters: primary mass A/p, mass-ratio 
q, orbital separation a, eccentricity e and metallicit y Z. Primary 
masses Mp are distributed assuming a power law JKroupa et al.l 
Il993h for the initial mass function (IMF). This is constrained by 
the o bservation of t he loc al l uminosity func tion and stellar den- 
sity of lWielen et"ai] h983h and |Poppeil3l98(]|) . A flat distribution is 
adopted for mass-ratio (0 < g < 1) and eccentricity (0 < e < 1), 
since these are not well constrained by observation. The distribu- 
tion of orbital separation a is assumed to be constant in logarithm 
for w ide binaries and falls off smoothly at close separations (Han 
Il998f) . We have assumed Z = 0.02 th roughout. All other input 
parameters are as given bv lYu & Jeffervl J2010,) . 

The present population synthesis simulation started with a 
sample of 10^ primordial MS+MS binaries, with total mass ~ 
1.05 X 10^ Mq, yielding a sample total of 4.93 x 10* DDs over 
all orbital periods of less than 100 yr during a sampling period of 
15 Gyr. 

In order to evaluate the computation of our DD sample, Fig.[T] 
shows (left) the orbital periods (log Porb/day) of each new-bom 
DD against primary mass (A/p) in the progenitors and (right) the 
time tivis from binary-star formation to the birth of its descendant 
DD. We define DDs in this sample as unevolved DDs. This figure 
shows the relation of new-born DDs and their progenitors. 



2.4 The contribution functions in the present simulation 

In this section, we investigate how the present number of DDs and 
the number of new-bom DDs depends on the star formation history 
of the thin disc by examining the contribution functions. We take 
the current age of the thin disc to be 10 Gyr. Star formation time 
tgf runs from to 10 Gyr. Figures|2] and [3] illustrate the variation 
of C (rate contributions to the present birth and merger rates) and 
C (number contributions to the present numbers and total merged 
numbers) with f^f . Equivalently, the figures show us the age distri- 
bution of present new bom, alive and merged DDs. Since the age of 



Figure 3. Same as Fig.|2]but for the contribution of star formation tgf to the 
present number (top panels) and the total merged number (bottom panels) 
of DDs, C. See ^^ior the details. 
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Figure 4. Birth rates u(t) (left panel) and merger rates (^(t) (left panel) 
and numbers n{t) (right panel) of different types of DD by integral of the 
corresponding birth and merger rates versus the age of the thin disc tdisc ■ 
u{t) and C^(t) for each type of DD is shown in the same colour, but C^(t) < 
u{t). Black points are for the nett total birth rate (i^(t) — C(*)) ^nd the total 
existing number of the DDs. 



a DD is defined as idisc — isf ( q2.3b . the age distribution is obtained 
by transforming the x-scale to idisc — isf • 

The top panels in Fig.[2] indicates that the present new-born 
He-l-He and CO-l-He DDs come mainly from early star formation 
tsf < 8 Gyr. This is not the case for CO-l-CO and ONeMg-l-X 
DDs. The top panels in Fig.|2]also show that the present new-born 
CO-l-CO and ONeMg-l-X come almost entirely from recent star for- 
mation; 8 < tsf < 10 Gyr 

This remarkable discrepancy is basically a result of iMS- Fig- 
ure[T] (right panel) illustrates that a massive MS-l-MS binary gen- 
erally takes less than a few Myr and 1 Gyr to form ONeMg-l-X 
and CO-l-CO DDs respectively, leading to their formation much 
more quickly than He-l-He and CO-l-He DDs. This also results in 
the birth rates of ONeMg-l-X and CO-fCO increasing dramatically 
at an early stage of disc evolution. Only after some time do CO-l-He 
and He-l-He DDs emerge and do their birth rates grow to reach the 
present values shown in the left panel of Fig.|4] Due to the quasi- 
exponential SF rate, the birth rates of ONeMg-l-X and CO-l-CO DDs 
start to decline after reaching a peak value. The long tMS for single 
He white dwarfs also results in a small contribution of star for- 
mation from 0—7.5 Gyr to the number of new-born ONeMg-l-He 
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DDs at the present. During the same period, there are no new-bom 
CO+CO DDs. 

Figure[3] shows that the contribution of star formation to the 
present numbers (top panels) and total merged numbers (bottom 
panels) of all types of DDs decreases monotonically as a function of 
tsf . This is because most DDs from all epochs survive to the present 
time due to their wide orbital separations. The number of He+He 
DDs from each epoch decreases with t^t more sharply than for 
CO+CO DDs, although early star formation provides more He+He 
DDs than CO+CO DDs. A similar situation arises for CO+He and 
ONeMg+He DDs. This result is consistent with stellar evolution 
and the assumed SF rate. 

The significance of computing the present number of DDs us- 
ing Eq.[6l which represents the sum of DDs arising from different 
star-formation epochs, is that it demonstrates the link between the 
SF history of the galaxy (or, at least, the thin disc in the present 
investigation) and the distribution of the properties of present-day 
DDs, which can be deduced from, for example, their gravitational 
wave signal. 

Figure|4] shows the variation of u, C^ and Wdd (Eq.[6ll of dif- 
ferent types of DD with age idisc- The individual properties of the 
curren t ndd DDs will be used to calculate the gravitational wave 
signal JYu&Jeffervll201ll) . 



Table 2. The contributions of different .star formation stages to the numbers 
of DDs at present day in the quasi-exponential star formation model in the 
thin disc. 



2.5 The structure of the thin disc and the local density of DDs 

The use of a realistic disc model is important in order to describe 
the distance distribution of white dwarf binary systems from the 
Sun. Sackett ( 1997) proposed a double exponential distribution. 
IPhleps et alj ((2000) derived three functions for the star density dis- 
tribution in their model of a thin disc plus thick disc (exponential 
+ exponential, hyperbolic secant + exponential, and squared hy- 
perbolic secant + exponential, respectively) from fits to deep star 
counts carried out in the Calar Alto Deep Imaging Survey. 

We here model the thin disc in the Galaxy using a squared 
hyperbolic secant plus exponential distribution expressed as: 



Pd{R,z) 



Mtn 



"^/''«sech^(-2//i,) MQpc"^ (10) 



where R and z are the natural cylindrical coordinates of the ax- 
isymmetric disc, ha = 2.5 kpc is the scale length of the disc, and 
hz — 0.352 kpc is the scale height of the thin disc. Mtn is the mass 
of the thin disc, which is determined by the star formation rate. We 
adopt the position o f the Sun to be -Rsun ~ 8.5 kpc, Zsun = 16.5 pc 
dFreudenreichll 19980 . We neglect the age and mass dependence of 
th e scale height . This thin- disc model is consistent with the model 
of iKlvpin et al.l ( l2002n and lRobin et al.l (|2003), and also in agree- 
ment with Hipparcos results and the observed rotation curve. 

From the SF rate, the total mass of stars in the thin disc at 
age 10 Gyr is Mtn « 5.2x10^° Mq. Combining the thin disc 
model and the mass of stars in the thin disc, the stellar den- 
sity in the solar neighbourhood is 6.27 x 10~^Mqpc~'^ for the 
thin disc, These values are co nsistent with the Hipparcos result. 



(7.6 ± 1.5) X 10"^ Mqpc "^ JCreze et al ."l998') and the dynami- 
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cal structure of the thin disc JKlvpin et alj 2002: Robin et al. 2003). 
The local density of DDs in the model is 1.1 



X 10"'' pc-^ 



3 IMPACT OF STAR FORMATION ON THE DD 
POPULATION 

The presence of unevolved DDs in Fig.[T] indicates the existence 
of a critical time when the first DD of each type was just born. 
For ONeMg+X, CO+CO, CO+He, and He+He DDs respectively, 
the times are 25 Myr, 50 Myr, 560 Myr, and 650 Myr. Figure|5]il- 
lustrates the contribution of different epochs of star formation to 
the present-day distribution of total mass and orbital periods of 
the DDs. We here distinguish the star formation for the current 
thin disc in six stages, which are (a) Sj isf/Gyr < 6 , (b) 
6 < ist/Gyr < 8 , (c) 8 ^ isf/Gyr < 9, (d) 9 < i,f/Gyr < 9.4, 
(e) 9.4 < isf/Gyr < 9.95, and (f) 9.95 ^ tsf/Gyr < 9.75, with 
the current disc age assumed to be 10 Gyr. We do not obtain any 
DD for star formation taking place after 9.975 Gyr. Figure|5] shows 
that the formation of very close compact binaries (log / > —2.5) 
is sensitive to star formation between 8 and 9.95 Gyr after the thin 
disc formed, which means that these DDs are most likely to be 
young. Their MS+MS progenitors formed between 50 Myr and 
2000 Myr ago. 

The total stellar mass formed during the time represented by 
each panel of Fig.[5]is 3.62, 0.85, 0.38, 0.15, 0.19, and 0.0087 x 
lO'^" Mq (a to f). The current number of DDs in the thin disc de- 
rived from each star formation epoch in the figure is given in ta- 
ble [2] These numbers indicate that, for current He+He and CO+He 
DDs, a large number (93.8% and 87.8%) have ages greater than 4 
Gyr, while only 71.9% and 67.8% of CO+CO and ONeMg+X have 
ages in the same range. A significant number (3.1% and 6.7%) of 
CO+CO and ONeMg+X DDs have been produced by the last 1 Gyr 
of star forination. The number of He+He and CO+He DDs from 
this period is negligible. However, IYu & Jeffervl ( 120111) show that 
DDs with ages less than 2 Gyr would contribute substantially to 
the amplitude of the gravitational wave signal in several frequency 
bands. 

According to the classification of the star formation stages and 
the critical time for the birth of DDs, we can see from Fig.[5]that the 
MS+MS progenitors of CO+CO DDs are formed before t^f ~ 9.95 
Gyr. The youngest CO+He DD has an age of about 560 Myr, but 
the majority of their MS progenitors formed >600 Myr ago. These 
results are consistent with stellar evolution calculations. 

Note that in the stellar evolution model a fraction of ONeMg 
white dwarfs become neutron stars and stellar-mass black holes 
due to accretion-induced collapse. These do not form type la su- 
pemovae and are not otherwise considered in our results. 



6 S.Yu& C.S. Jejfery 




1-2-3 4 3 2 10-1-2-3 

3.0r 



2.5 
2.0 
1.5 
1.0 
0.5 







1" 


-1- 


He+He WDs 




X 


CO+Hc WDs 




A 


CO+CO WDs 


SFat 


O 


ONe WDs 


9.95-9.975 Gyr 



S^Jdays) 



og P^^ (days) 



4 3 2 10-1-2-3 

log P ,, (days) 



Figure 5. The influence of different stages of star formation on the orbital periods and the total mass (primary mass -I- secondary mass) of present DDs. 
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Figure 6. The rate contribution function of star formation epochs to the 
present birth rates (C, .new , top) and merger rates (C7* ,mer , bottom) of DDs 
with respect to star formation time tgf in sample. 



4 COMPARISON OF DIFFERENT STAR FORMATION 
MODELS 

We have simulated the present DD population with a quasi- 
exponential declining SF rate. In order to see the influence on the 
DDs from different SF models, we have arbitrarily constructed an- 
other three simple SF models and calculated the birth rates, merger 
rates, and the number of DDs. In addition to the quasi-exponential 
SF model described above, these are : 
Instantaneous SF. A single star burst takes place only at the for- 
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the different star formation models. 




,/Gyr 
Figure 9. As Fig.[8]for CO+He DDs. 



10"' 


/^^ : 


10-' 


" f.^?x^'^' 


10-= 


" / /' 1 


- 


10-' 


i' 





0.1 1 10 

^a:./Gyr 

Figure 10. As Fig.Hfor CO+CO DDs. 





^10-' 



10- 




instantaneous 

constant 

enhanced 

quasi-exp 



0.1 



'.. /Gyr 



10 



Figure 11. As Fig.[8]for ONeMg+X DDs. 



We take to = 391 Myr, so the resulting SF rate in this model is 
132.9 Mq yr"'^ on average from ist = to 391 Myr. There is no 
subsequent star formation. 
Constant SF. A constant SF rate at an average rate of 



S = 5.2 M0 yr"\ sS tdisc < 10 Gyr. 



(12) 



Enhanced SF. An enhanced SF rate during tdisc from 8 to 10 Gyr, 



f 9, < tdisc < 1 Gyr, 
S = <^ 4.4, 1 < idisc ^ 8 Gyr, 
I 6, 8 < tdisc ^ 10 Gyr. 



(13) 



Each model produces a present stellar mass of 5.2 x 10 Mq 
in the thin disc with present age f disc = 10 Gyr in ag reement with 
the dynamic mass of the thin disc JKlvpin et al.ll2002l) . 

Figures[6] and [7] show the star-formation contribution func- 
tion C, in our sample to the birth rates and merger rates of the 
DDs and the SF rate for each SF model. These two variables are 
needed to calculate the birth rates and merger rates in the thin disc 
in Eq.|7] The many small oscillations on the star-formation contri- 
bution functions in Fig.|6](and Fig.|2j are due to the statistical noise 
of the Monte Carlo simulations. 

Figures [8l to 1 1 1 1 show the birth rates (i/) and merger rates (Q 
of the four types of DD in the four models. We see that the instan- 
taneous SF differs greatly from the other three models. As a result 
of the impact of a high SF rate at the formation of the thin disc, f 
and (^ in the instantaneous SF model increase very sharply to reach 
a maximum value several times greater than in all other models, f 
and C subsequently decrease slowly for the He-l-He DDs since they 
have a long iMS. In the cases of the other three types of DD in the 
instantaneous SF model, the decline of v and (^ is faster than He-l-He 
DDs due to the shorter tMS. When tdisc > 2.3 Gyr, roughly the 
longest tMS of CO-fCO DD, there are no more new-bom CO-l-CO 
DDs. ONeMg-l-X DDs show a reduction in i^ and ^, except for the 
ONeMg-l-He DDs which make a residual contribution to u as the 
disc continue to evolve. Again, some oscillations are seen in the 
birth and merger rates in the instantaneous SF model. These are 
due to statistical noise from the contribution functions produced by 
the Monte Carlo simulations, and numerical noise from the quadra- 
ture of Eq.|7] 

A big difference for u and (^ in the other three models can 
also be seen from the figures. Enhanced SF causes the present u 
and C of CO-l-CO and ONeMg-l-X DDs to be slightly higher than 
constant SF, while the quasi-exponential SF can produce a higher 
u and C of CO-fCO and ONeMg-l-X DDs than enhanced SF when 
tdisc < 4 Gyr. When tdisc > 4 Gyr, we find f and ( of CO-l-CO and 
ONeMg-l-X DDs decrease dramatically in the quasi-exponential SF 
model, while the values of ly and ^ in the constant SF model remain 
almost constant and higher than the v and ^ in the quasi-exponential 
SF model, up to tdisc = 10 Gyr. The significant decreases and 
increases of the values of i/ and ^ of CO-l-CO and ONeMg-l-X DDs 
in the enhanced SF model are the response of these variables to 
sudden changes in the SF rate. The present values of i/ and ( are 
similar because we adopt a similar value of the average SF rate at 
the present age of the thin disc in these three models. 

The present number and the total merger number of different 
types of DDs are shown in Figs. [12] to [15] A similar evolutionary 
history of the numbers for various types of DD is generated by the 
continuous SF models since the average values of the SF rate are 
similar in these SF models. The instantaneous SF model produces 
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Table 3. As table[2]except for the instantaneous star fomiation model. 
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Figure 13. As Fig.[T2]for CO+He DDs. 



a slightly higher present numbers of DDs than the other three SF 
models, and has a distinctly different evolution. 

A delay time, which represents the time from when the first 
DD was born to when the first DD merged, can certainly be found 
in each model if we compare the left panel and the right panel in 
each figure. The delay time can not be less than the shortest life 
time toD of a DD, and is affected significantly by stellar evolution 
models and population synthesis parameters (e.g. the IMF). Since 
we use the same stellar evolution model and population synthesis 
parameters for each single epoch of the SF in our simulations, the 
delay time is a constant for each SF model. We list the delay time 
for different types of DD in Table|6] 



0-0.391 Gyr 
0.391-10 Gyr 



He+He 

112790783 





CO+He 

35296981 




CO+CO 

74199554 




ONeMg+X 

7979607 





0-10 Gyr 



112790783 35296981 74199554 



7979607 



5 SUPERNOVAE 

We have been concerned about the details of birth rate and number 
of DDs in the present simulation. However, we suggest that Eq.|7] 
(to calculate the birth rate) and Eq.|6]and[8](to calculate the number) 
are valid for almost any type of star. The only thing required is to 
determine the corresponding contribution of each SF epoch to the 
numbers of the stars in question with respect to the variation of the 
SF {i.e. tsi). In terms of the determination of the contribution, it is 
important to determine all possible formation channels. Otherwise 
we should specify that the birth rate or number correspond to a 
specific formation channel. We would like to take two important 
examples, type la supernovae (SNIa) and core collapse supernovae 

(ccSN, including type Ib/Ic and type UX 

Previous studies (e. g. Ibe n & Tutukovl (ll984h:lTutukov et al.l 
( 1992); Iben & Tutuko^ i ll997h : lYungelson& Livid l ll998ll20o'ol) : 
IDahlen et aljt2004i) ) indicate that mergers of DDs with a total mass 
exceeding the Chandrasekhar mass limit (Mch) would be one for- 
mation channel to produce SNIa. We here investigate the birth rate 
of SNIa from the mergers of DDs, and we constrain the progenitors 
to be those DDs with total mass > Afch ~ 1.38 Mq which can 
merge in the time tdisc to tdisc + <5tdisc. 
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Table 4. As table|2]except for the constant star formation model. 



Figure 14. As Fig.[T2]for CO+CO DDs. 
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Table 5. As table[2]except for the enhanced star formation model. 
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Table 6. The delay time, which is from the time the first DD was born to 
the time the first merged DD appears, for each type of DD. 



He+He CO+He CO+CO ONeMg+X 



delay time (Myr) 724.03 23.94 



28.65 



22.39 



The birth rates of SNIa for the four SF models are shown in 
the left panel in Fig.[T6l and the rate contribution function from 
each single epoch in the four models is shown in the right panel in 
the same figure. The calculation of the rate contribution function is 
from the number of DDs satisfying the conditions to become SNIa 
in our simulation. Note that the curve in the right panel is also the 
age distribution of SNIa, taking into account the SF rate in the thin 
disc. This figure indicates that both old and young DDs contribute 
to the present SNIa rates. 

From the left panel, we see that the birth rate of SNIa from the 
instantaneous SF model reaches its maximum value very quickly 
with the evolution of the thin disc. The maximum value is signifi- 
cantly higher than that in another three SF models since the initial 
SF rate in the instantaneous SF model is more than ~30 times that 
in other SF models. The enhanced SF model produces the high- 
est present rate of SNIa in the four SF models. This is because the 
present SF contributes the majority of mergers. Again, from Fig.ll6l 
continuous SF is also the reason why the rate of SNIa keeps grow- 
ing with the the evolution of the thin disc. The present-day SNIa 
rate is about 2.19, 4.92, 5.07, and 4.34 xlO"** yr"^ for the in- 
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Figure 16. The rates of type la supemovae from double degenerate mergers 
channel in different star formation models (left panel), and the rate contri- 
bution function to the prensent-day (tdisc = lOGyr) supemovae rates in 
the thin disc (right panel). 



stantaneous, constant, enhanced, and quasi-exponential SF models 
respectively. 

We have to emphasize that we here only take into account 
the merger channel for the rate of thin-disc SNIa. This leads to 
our results being less than the occurrence rate inferred from ob- 
servations, ~ 4 X 10~^ yr~^ JCappellaro et aI.lll997D , as a few 
other form ation channels would contribute to the observed occur- 
rence rate ( Yungelso n & Lividl200d : iHan & Podsiadlowskill2004l : 
iHachisu et al. 2008) . 

We note that the occurrence (or birth) rate of SNIa at tdisc = 
10 Gyr in our simul ations is sma ller than in som e other recent the - 
oretical studies (e.g. iRuiter et all (2009); Menne kens et al.l ( |2010|) ) 
by a factor of 1.5 — 4. Since all these studies (including ours) adopt 
the same IMF for the primary main sequence stars, the same flat 
distribution for initial mass ratio, and the same initial eccentricity 
distribution, we infer that the differences between our study and 
others are mainly caused by the initial distribution of orbital sepa- 
rations and by the treatment of mass transfer and mass loss during 
binary ev olution. 

Both lRuiter et al.l ( |2009|) and lMennekens et al.l(l2010l) a dopt a 
distribution of the orbital separation ^ oc a~^ (Abt 1983) when 
a ^10^ R p,, while we use a distribution of oc a"'^ ( Griffin 19851 ; 
lHai]| 19981) for a ^ 10 Rq, and oc a"^ for 10 :? a sj 5.75 x 10*^ 
7?0. This may lead to the generation of more close DDs in their 
models than our model. 

For co mmon-envelope ejection, we adopt the 7-algorithm 
JNelemans & Tout.2005.) a nd fix 7 = 1.5. Other studies use the q- 
algorithm JWebbinkll 19840 . Increasing the common-envelope ejec- 
tion efficiency (for example, inreasing qA = 0.5 — >■ 1.0) results in 
more binaries going through a double common-envelope phase to 
become close DDs, and hence forming potential SNIa progenitors. 
A close comparison between the studies implies that the common- 
envelope ejection efficiency is lo wer in our calcula tion than in oth- 
ers. In the most extreme case, the lYungelsonl J20I (J) SNIa birth rate 
from the DD merger channel is higher than ours by a factor of '--^lO 
because they adopted a very high common-envelope ejection effi- 
ciency (aA = 2.0) and a higher SF rate (~8 Mq^x^^'). 

A significant difference between the studies concerns the 
star-formation his tory for early-type (e lliptical) galaxies. Both 
iRuiter et al.l j2009l) and lMennekens et al.l 12010) approximate "in- 
stantaneous star formation" by a (5-function star-burst. Our approx- 
imation of "instantaneous SF" assumes that the initial star-burst 
lasted a few million years at a constant rate. This leads to the peak 
value of the SN Ia rate occur ring rather later i n our models tha n 
in the studies in IRuiter etakl 12009) and Mennekens et alj ( l2010h . 
However, none of these approximations may reflect the true SF 
history of early-type galaxies; recent studies indicate th at bright 
early- type galaxies show signs of current star formation l IYi et al.l 
1 20051) and, by inferrence, may have experienced intermittent star- 
burst episodes. 

The case of ccSN is easier than that of SNIa. Their progeni- 
tors, main sequence stars with mass approximately greater than 8 
Mq, evolve very quickly on a timescale <20 Myr. This is even 
shorter than the tMS of ONeMg+X (see Fig.[T]l, so all the new- 
born ccSN results from SF after a few tens Myr (see criterion in 
i]2.3l ). Then the SF rate can be regarded as a constant ^(tdisc). 
From Eq.|7] we have the birth rate of the supemovae SNR{tdisc) ~ 
[A''ccSN/mtotai] • ^(tdisc), where iVccSN is the number of stars go- 
ing supernova in a calculation of a sample of main sequence stars, 
w^totai represents the total mass of the sample , and tdisc = 10 Gyr. 
Using the same parameters as in the simulation in this paper, we 
obtain iVccSN/mtotai ~ 4.5 x 10"^ Mg^ and so SNrIw) ^ 
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l.JiA, 2.7 and 1.58 century ^ for the three continuous SF mod- 
els. These va lues are consisten t with the observations of 1.2 — 
3 century"^ jS mith et al."l978'; 'van den Bersh & McClure"l994'; 
iMcKee & Wmi ams 1997; Timmes et al. 1997; Diehl et al. 2006) in 
terms of the mass of ^'^ Al, the number of massive stars in the local- 
ized HII region and the extragalactic supernova scaled calculation. 



6 DISCUSSION 

Our simulation is based only on theory with some simple assump- 
tions. The final outputs, the rates and numbers, are derived from 
two key inputs which are the SF rate and the contribution func- 
tions for the evolved stars. On the other hand, if we know the rates 
and numbers of objects in a galaxy, especially some types of exotic 
stars, we are able to do a deconvolution by setting the observed rate 
as an input parameter to derive the SF rate at specific epochs. In 
fact, this method is used to obtain the current SF rate via the ob- 
servation of supemovae JHeger et alj|2003l ; lOiehl et al.ll2006r) and 
a global SF rate by obse rvations of late type main sequence stars 
JRocha-Pinto et alj|2000h . The contribution functions provide an 
important link between SF history and the diverse stellar compo- 
nents in different galaxies, since they reflect the sensitivity of a 
group of the same type of stars to the SF history of a galaxy. 

The time grid Stst (or Stdisc) can cause some uncertainties. 
For the quasi-exponential declining SF rate, another computation 
with a low-resolution time grid of St^t = 0.026 results in posi- 
tive deviations from the high resolution time grid (Stst — 0.00867) 
of 1.03%, 0.82%, 0.04% and 0.21% for the birth rates of He-l-He, 
CO+He, CO-fCO and ONeMg-l-X DDs respectively; the deviations 
of the numbers are 3.97%, 2.53%, 1.13% and 0.87% for the differ- 
ent types of DDs. In the new grid, the deviations of the numbers 
calculated from Eq.|6] compared with Eq.[8]are 4.7%, 2.8%, 1.5% 
and 1.4% for each type of DD. The time error for different stellar 
evolution phase in our simulation is less than 4.2%. 



7 CONCLUSION 

In this paper we have investigated the relation between the evolu- 
tion of birth rates, merger rates and numbers of different types of 
DD and the SF history in the thin disc of the Milky Way Galaxy. 
By analysing how a quasi-exponentially declining SF rate influ- 
ences the rates and observed numbers of DDs with respect to the 
evolution of the disc and SF, we find that SF between and ~8 
Gyr dominates the present rates and total numbers of He+He DDs 
and CO+He DDs. Similarly, the current numbers of CO+CO and 
ONeMg-l-X DDs mainly come from early SF, although the later SF 
(e.g. 4 to 8 Gy) contributes more CO+CO than it does for He+He 
which are the two largest DD population in the thin disc. 

However, the present birth and merger rates of CO+CO and 
ONeMg+X DDs are strongly governed by the recent star formation 
(i.e. 8 to 9.95 Gyr). SF before ~7.5 Gyr does not contribute to 
the present birth rates of CO+CO ang ONeMg+X DDs, but it has 
some contribution to the merger rates. More importantly, with the 
SF history related population synthesis approach in this paper, not 
only are we able to determine the rates and numbers of DDs, but 
we can also obtain the distributions of properties of current DD 
population from different stages of SF. 

We have compared the impact of different SF models, namely 
the instantaneous, the constant, the enhanced, and the quasi- 
exponential SF, on the rates and numbers of DDs and the rates of 



two types of supernovae. The evolution tracks of the rates and num- 
bers from the four models are quite different. A distinct difference 
can be found between the instantaneous model and the other three 
continuous SF models. This model gives historical rates and num- 
bers obviously higher than the other three models. However, the 
present rates of DDs from this model are apparently lower than the 
other three models. 

In addition to the DDs, we have calculated the rates of SNIa 
and ccSN. The evolution of the rates of SNIa is basically similar 
in the four SF models, but the instantaneous model can produce 
a higher rate in the past because of the very high SF rate at the 
formation of the thin disc. The present rates of SNIa are 2.19, 4.92, 
5.07, and 4.34 x 10"'* yr~^ for the instantaneous, the constant, the 
enhanced, and the quasi-exponential SF models respectively. The 
rates of ccSN from all four SF models, ~1.5 to 3 century"^, are 
consistent with the observations. 



ACKNOWLEDGMENTS 

The Armagh Observatory is supported by a grant from the Northern 
Ireland Dept. of Culture Arts and Leisure. The authors are grateful 
to Prof. Lev Yungelson for useful suggestions, and to colleagues at 
the Armagh Observatory for their stimulating discussion. We thank 
the referee for the useful suggestions and comments. 



REFERENCES 

Abt H. A., 1983, ARA&A, 21, 343 

Cappellaro E., Turatto M., Tsvetkov D. Y., Bartunov O. S., Pollas 
C, Evans R., Hamuy M., 1997, A&A, 322, 431 

Creze M., Chereul E., Bienayme O., Pichon C, 1998, A&A, 329, 
920 

Cropper M., Harrop-Allin M. K., Mason K. O., Mittaz I. P. D., 
Potter S. B., Ramsay G., 1998, MNRAS, 293, L57 

Dahlen T., Strolger L., Riess A. G., Mobasher B., Chary R., Con- 
selice C. J., Ferguson H. C, Fruchter A. S., Giavalisco M., Livio 
M., Madau P, Panagia N., Tonry J. L., 2004, ApJ, 613, 189 

Diehl R., Halloin H., Kretschmer K., Lichti G. G., Schonfelder 
v.. Strong A. W., von Kienlin A., Wang W., Jean P., Knodlseder 
J., Roques J., Weidenspointner G., Schanne S., Hartmann D. H., 
Winkler C, Wunderer C, 2006, Nature, 439, 45 

Evans C. R., Iben I. J., Smarr L., 1987, ApJ, 323, 129 

Farmer A., Roelofs G., 2010, ArXiv e-prints 

Freudenreich H. T, 1998, ApJ, 492, 495 

Gokhale V., Peng X. M., Frank J., 2007, ApJ, 655, 1010 

Griffin R. F, 1985, in P P Eggleton & J. E. Pringle ed., NATO 
ASIC Proc. 150: Interacting Binaries The distributions of periods 
and amplitudes of late-type spectroscopic binaries, pp 1-12 

Hachisu I., Kato M., Nomoto K., 2008, ApJ, 683, L127 

Han Z., 1998, MNRAS, 296, 1019 

Han Z., Podsiadlowski P, 2004, MNRAS, 350, 1301 

Heger A., Fryer C. L., Woosley S. E., Langer N., Hartmann D. H., 
2003, ApJ, 591, 288 

Hurley J. R., Pols O. R., Tout C. A., 2000, MNRAS, 315, 543 

Hurley J. R., Tout C. A., Pols O. R., 2002, MNRAS, 329, 897 

Iben Jr. I., Tutukov A. V., 1984, ApJ, 54, 335 

Iben Jr. I., Tutukov A. V., 1997, ApJ, 491, 303 

Klypin A., Zhao H., Somerville R. S., 2002, ApJ, 573, 597 

Kroupa P, Tout C. A., Gilmore G., 1993, MNRAS, 262, 545 

Majewski S. R., 1993, ARA&A, 31, 575 



Star formation history and double degenerates 1 1 



Marsh T. R., Nelemans G., Steeghs D., 2004, MNRAS, 350, 113 

McKee C. R, 1989, ApJ, 345, 782 

McKee C. R, Williams J. R, 1997, ApJ, 476, 144 

Mennekens N., Vanbeveren D., De Greve J. R, De Donder E., 
2010, A&A, 515, A89+ 

Nelemans G., Tout C. A., 2005, MNRAS, 356, 753 

Nelemans G., Yungelson L. R., Portegies Zwart S. R, 2001, A&A, 
375, 890 

Nelemans G., Yungelson L. R., Portegies Zwart S. R, 2004, MN- 
RAS, 349, 181 

Phleps S., Meisenheimer K., Puchs B., Wolf C, 2000, A&A, 356, 
108 

Popper D. M., 1980, ARA&A, 18, 115 

Ramsay G., Hakala P., Wu K., Cropper M., Mason K. O., Cordova 
R A., Priedhorsky W., 2005, MNRAS, 357, 49 

Robin A. C, Reyle C, Derriere S., Picaud S., 2003, A&A, 409, 
523 

Rocha-Pinto H. J., Scalo J., Maciel W. J., Flynn C, 2000, A&A, 
358, 869 

Ruiter A. J., Belczynski K., Benacquista M., Larson S. L., 
Williams G., 2010, ApJ, 717, 1006 

Ruiter A. J., Belczynski K., Fryer C, 2009, ApJ, 699, 2026 

Sackett P D., 1997, ApJ, 483, 103 

Smith L. R, Biermann P, Mezger P G., 1978, A&A, 66, 65 

Timmes R X., Diehl R., Hartmann D. H., 1997, ApJ, 479, 760 

Tutukov A. v., Yungelson L. R., Iben Jr. I., 1992, ApJ, 386, 197 

van den Bergh S., McClure R. D., 1994, ApJ, 425, 205 

Webbink R. R, 1984, ApJ, 277, 355 

Wielen R., JahreiB H., Kruger R., 1983, in A. G. D. Philip & 
A. R. Upgren ed., lAU CoUoq. 76: Nearby Stars and the Stel- 
lar Luminosity Runction The Determination of the Luminosity 
Function of Nearby Stars, pp 163-170 

Yi S. K., Yoon S.-J., Kaviraj S. e. a., 2005, ApJ, 619, LI 11 

Yu S., Jeffery C. S., 2010, A&A, 521, A85-^ 

Yu S., Jeffery C. S., 2011, Submitted 

Yungelson L., Livio M., 1998, ApJ, 497, 168 

Yungelson L. R., 2010, Astronomy Letters, 36, 780 

Yungelson L. R., Livio M., 2000, ApJ, 528, 108 



